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Abstract 

The spin-dependent twist-3 and 4 nucleon matrix elements d^ and describing 
power corrections to the Bjorken and Ellis-Jaffe sum rules are computed in the 
instanton vacuum. A systematic expansion in the small packing fraction of the 
instanton medium, p/R <C 1, is performed. We find that the twist-3 matrix element 
d^ is suppressed [of order (p/R)^], while the twist-4 matrix element is of 
order unity. Numerically, d^ <C f( 2 \ The small value of d^ ~ 10~ 3 obtained from 
instantons is consistent with the recent E143 measurements of the structure function 
g2, where d^ enters at the same level as the twist-2 contribution. 
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In QCD, the quantitative description of deep-inelastic scattering and other related 
processes comes down to a study of matrix elements of certain composite operators in the 
nucleon state. The well-known parton distributions, which parametrize the non-power 
suppressed part of the structure functions, are given by matrix elements of operators of 
twist 2. In addition, there is a large number of observables related to operators of non- 
leading twist (> 2). One example are power corrections to the Bjorken sum rule |], @. To 
order l/Q 2 : 

£dxgr(x,Q>) = iog + ^(og + 44g + 4/$), (1) 

where a^\a^ denote the matrix elements of twist-2 operators, with a^l = (<7a)ns 5 while 
and parametrize, respectively, the proton matrix element of the twist-3 spin-2 
and the twist-4 spin-1 operators: 

(PS\0$* + Og?\PS) - traces 
= 2M N djg [2P a P p S^ - P^P p S a - P a P' 1 S f3 + {a^ (3)- traces] , (2) 

(PS\Of s a \PS) = 2M% r fj$St>, (3) 

0&{x) = g^(x)r 3 rG^(x)^(x). (4) 

Knowledge of the values of S 2 ^ and would allow one to extend the QCD analysis of 
the Bjorken and Ellis- Jaffe sum rules to low values of Q 2 , where data with good statistics 
is available, which could be of interest e.g. for an accurate determination of as- Similar 
expressions can be written for the flavor-singlet part of g\ (Ellis-Jaffe SR), as well as for 
the GLS sum rule for F 3 , see Refs.JI], |3|. 

It should be stressed that the parameters S 2 \ f^ 2 \ and the like, are universal char- 
acteristics of the nucleon, whose significance extends far beyond power corrections to g±. 
For instance, the twist-3 matrix element S- 2 ^ 1 appears at leading twist level (i.e. not power 
suppressed relative to twist-2) in the third moment of the structure function g 2 M' 

jf dxx 2 9 r n (x,Q 2 ) = -iog + + oQ_) (5) 

Also, the twist-4 matrix element can be related, via the heavy-quark expansion, to the 
charm contribution to g\ 0. Thus, these matrix elements can in principle be determined 
by combining results of different measurements. In order to compute them from first 
principles, one needs a theory of the nonperturbative effects giving rise to quark-gluon 
correlations in the nucleon. 

There is by now overwhelming evidence — both from phenomenology as well as from 
lattice simulations || - for an important role of instanton-type vacuum fluctuations in 
determining the properties of the low-energy hadronic world ||. The so-called instanton 
vacuum approximates the ground state of Yang-Mills theory by a medium of instantons 
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and antiinstantons, with quantum fluctuations about them, which is stabilized by instan- 
ton interactions J7|. The coupling constant is fixed at a scale of the order of the inverse 
average size of instantons in the medium, ~ 600 MeV. The most important property 
of this picture is the diluteness of the instanton medium, meaning that the ratio of the 
average size of instantons in the medium to their average distance, R, is small: p/R~ 1/3 
|6|, [7j. The existence of this small parameter makes possible a systematic analysis of 
non-perturbative effects in this approach. 

In particular, the instanton vacuum explains the dynamical breaking of chiral symmetry 
II |§ . It happens due to the derealization of the fermionic zero modes associated with the 
individual (anti-) instantons in the medium. Using the l/iV c -expansion one derives from 
the instanton vacuum an effective low-energy theory, whose degrees of freedom are pions 
(Goldstone bosons) and massive "constituent" quarks. It is described by the effective 
action [HI 



S cS = J d A x ip(x) - MF(9)e !75T ° lOW F(i) 



1>{x). (6) 



Here, M is the dynamical quark mass generated by the spontaneous breaking of chiral 
symmetry; parametrically 

MP ~ (7) 

and F(k) is a form factor proportional to the wave function of the instanton zero mode, 
which drops to zero for momenta of order k ~ p~ l . Mesonic correlation functions, com- 
puted either with the effective action, Eq.@, using the l/iV c -expansion ||, or by more 
elaborate numerical simulations ||, show excellent agreement with phenomenology. The 
nucleon can be described, within the l/iV c -expansion, as a chiral soliton of the effective 
theory flO|| . This field-theoretic description of the nucleon gives a good account not only 



of its static properties, but also of the twist-2 parton distributions [11 . 

Given the general success of the instanton vacuum in describing hadronic properties, 
it is natural to apply this approach to the calculation of matrix elements of higher-twist 
operators such as Eq.(0) and (|3|). A general method for computing matrix elements of 



quark-gluon operators in the instanton vacuum has been developed in Ref. ||12||. Due to 
the diluteness of the instanton medium the dominant contribution to the matrix element 
is obtained by substituting for the gluon operator in Eqs.(|2|) and (|3|) the classical gluon 
field of a single (anti-) instanton. The interaction of the instanton with the fermions 
through the zero modes then allows one to replace the original QCD operator (normalized 
at p = p _1 ) by an effective operator, expressed in terms of the quark and pion fields 
of the effective theory. It was shown in Ref.[Q that identities following from the QCD 
equations of motion are fully preserved in this approach. For the operator Eq.(|j) the 
effective operator, after integration over the instanton coordinates, takes the form: 



■ l O n a/37 (x) = ip(x)T y7 a ^(a;) x ^ / d z fp^^x - z) 



X 



^{z)^a, v FCd)e^ a{z) FCd)m 
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where the function //? 7)AtI ,(x) describes the (anti-) self dual field of the (anti-) instanton in 
singular gauge, 



G%y( x )l(I) — ±G a n Jx)ui) — ~ffa,ya>{x) < , (9) 
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f/3 7 M x ) - _^ 2)2 {~^r 5 ^ + ~^ 5 ^ ~ 2 6 ^ u ) ' ^ 

In Eq.(||) the gluon field of the QCD operator has been replaced by a chirally-odd 
instanton-induced quark vertex proportional to the dynamical quark mass, M. Note 
that, as in the effective Lagrangian, Eq.(|]), chiral invariance is preserved due the presence 
of the pion field. 

The matrix element of the effective operator, Eq.(Q), can now be computed in the 
effective low-energy theory. At this point characteristic differences between the different 
spin projections of the operator emerge. A parametrically large contribution can only 
come from contractions of the four-fermionic operator, Eq.(^), with closed quark loops, 
in which the quark momenta run up to the ultraviolet cutoff, p~ l (see Ref.@ for details). 
One finds that in the case of the twist-4 matrix element, f^ 2 \ Eq.(^), the loop is oc pr 2 
( "quadratically divergent"), while in the case of the twist-3 matrix element, S- 2 \ Eq.(|J), it 
is only of order logMp, ("logarithmically divergent"). As a result, we obtain a parametric 
difference of the twist-4 and 3 matrix elements: 

^° 
R 

R 



twist 4: /( 2 ) ~ (Mp)° 
twist 3: ~ (M p) 2 log Mp ~ (J^j log 



This remarkable result can be traced back to two circumstances: the 0(4) symmetry of 
the field of a single (anti-) instanton, and the diluteness of the instanton medium. 

The parametric difference between the twist-3 and 4 matrix elements also reflects itself 

(2) 

in the numerical values. For the flavor non-singlet nucleon matrix element (the NS 
part is leading in the l/iV c -expansion) we obtain a value of —0.1, while <%s was estimated 
to be of the order of 10~ 3 .[] Thus, the instanton vacuum implies a clear numerical ordering 
of higher-twist matrix elements: \S 2 ^\ |/^|. Again we stress that this ordering is not 
accidental, but a direct consequence of qualitative properties of the instanton medium. 
This is at variance with the QCD sum rules estimates, and also of bag model calculations, 

(2) (2) 

which typically give values of and %g of the same order, see Table 1. (For a critical 
discussion of the QCD sum rule calculations see Ref.|i"7|.) Note also that the instanton 



(2) 

result for agrees in sign with the QCD sum rule results, but not with the bag model. 
It is interesting that the small value of dP> suggested by the instanton vacuum seems 



to be in agreement with recent measurements of #2 by the E143 collaboration. Ref . [jig] 



^^An accurate calculation of in this approach would require to take into account two-instanton 
contributions in the effective operator; this calculation is presently being done. The order-of-magnitude 
estimate in Ref.|| is based only on the one-instanton contribution. 
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Instanton vacuum 
QCD sum rules |13 
QCD sum rules [14 
Bag model 15] 
E143 16] 


1 


-1.0 
-2.0 
-0.7 
1.1 


~ 0.01 
0.7 
0.7 
0.6 

0.024 ±0.78 


0.4 
1 
1 

5 
5 



Table 1: Numerical results for the flavor-nonsinglet spin-dependent twist-4 and 3 matrix 
elements /^g = 3(/( 2 ) — ffl) and rf^s, Eqs.(||]) and (|3|). The instanton result for d!^\ is an 
order-of-magnitude estimate (see the text). 

reports a value of d^s = 3(c^ ) 2 ' ) — dffi) = (2.4±78) x 10~ 3 , which agrees well with the order- 
of-magnitude estimate of Ref.[Q (note that df> and individually are of comparable 
magnitude, see Ref.[16|]). This value is much smaller than our estimate for /^ g , which 



is of order 10 _1 . Unfortunately, the present data are not sufficient to unambiguously 
extract from power corrections with sufficient accuracy in order to make a quantitative 
comparison of and d^ ||18 |. It would be interesting to see if by imposing the condition 



|d (2) | < |/ (2) | from the start one could extract more accurately. 

To summarize, we have shown that the instanton vacuum with its inherent small param- 
eter, p/R, implies a parametrical (and numerical) hierarchy of the spin-dependent twist-3 
and 4 matrix elements: \d^\ <C |/^|- This qualitative prediction can be confronted with 
structure function data. Extending this analysis further (e.g. to higher moments and to 
other structure functions) one could use the rich information available from DIS to learn 
about properties of non-perturbative vacuum fluctuations. 

It is a pleasure to thank D.I. Diakonov, V.Yu. Petrov, P.V. Pobylitsa for many inter- 
esting discussions, and K. Goeke for encouragement and support. 
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